RPA CORRELATIONS AND NUCLEAR DENSITIES IN 
RELATIVISTIC MEAN FIELD APPROACH 



N. VAN GIAli'2, H.Z. LIANG^^^.a^ MENG^ 

^InsUtut de Physique Nudeaire, CNRS, UMR8608, F-91406 Orsay, France, E-mail: 

nguyenQipno. in2p3.fr 
Universite Paris-Sud, F-91406 Orsay, France 
^School of Physics, Peking University, 100871 Beijing, China 

' (Dated: February 5, 2008) 

o ■ 
o 

. Abstract. The relativistic mean field approach (RMF) is weU known for describing 

' accurately binding energies and nucleon distributions in atomic nuclei throughout the nuclear 

^ \ chart. The random phase approximation (RPA) built on top of the RMF is also a good 

framework for the study of nuclear excitations. Here, we examine the consequences of long 
fT^ ' range correlations brought about by the RPA on the neutron and proton densities as given 

by the RMF approach. 

^ Key words: Relativistic Mean Field, Random Phase Approximation, nucleonic densi- 

^ : ties. 

^ 1 INTRODUCTION 

o ■ 

■ The self-consistent mean field approach has become a powerful theory for 

understanding the structure of atomic nuclei at a microscopic level. The 
building blocks are structureless nucleons interacting via two-body effective 
interactions. At lowest order of the theory, correlations are neglected and the 
^ ■ nucleus is described as a system of A particles moving independently in a mean 

\ field created self-consistently by the two-body interactions. This mean field is 

determined by applying a variational principle to the total energy. 

The above picture is known as the static Hartree-Fock (HF) theory. A 
natural extension of it is the time-dependent Hartree-Fock (TDHF) approach, 
\ whose linearized version coincides with the well-known random phase approxi- 

mation (RPA)[lJ. In this way, long-range correlations of the particle-hole type 
are introduced. They will bring corrections to various properties predicted 
at the HF level, such as total energies, single-particle spectra and occupation 
probabilities, densities, etc... In this work, we want to focus on the effects 
of RPA long range correlations on the nuclear densities. This question has 
been studied several decades ago in the framework of HF-RPA calculations 
done with the Gogny effective interaction [2]. More recently Dupuis et al.[3j 
have re-examined the influence of RPA correlations on nuclear densities and 
potentials, using again the Gogny interaction. 
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However, besides the non-relativistic HF and RPA framework there is 
the relativistic counterpart represented by the relativistic mean field (RMF) 
approach [H El E] and the relativistic RPA (RRPA) built upon the RMF[7, 8j. 
The RMF approach consists in treating a relativistic meson-nucleon effective 
Lagrangian at the level of the Hartree approximation (no exchange) and no-sea 
approximation (i.e., the Dirac sea is considered empty). Effective Lagrangians 
with density-dependent meson-nucleon couplings have been determined in or- 
der to give an excellent description of nuclear ground states in RMF[9l|6] as 
well as collective excitations in RMF-RRPA[8j. It is thus worthwhile to ex- 
amine the influence of RRPA correlations on the ground state properties and 
to evaluate how much the RMF predictions will be affected. This is what we 
will do in the rest of this paper, concentrating the discussion on the nuclear 
densities. 



2 THEORY 

2.1 REMINDER OF RRPA 

The starting point is an effective Lagrangian where the nucleon-nucleon in- 
teraction is mediated by the exchange of effective mesons (scalar-isoscalar a, 
vector-isoscalar uj, vector-isovector p) and photons. The meson-nucleon cou- 
plings are assumed to be density dependent. In general, the parameters are 
adjusted so as to give a good description of nuclear ground states throughout 
the periodic chart in the RMF approach, i.e., in static Hartree plus no-sea 
approximation. Having solved the RMF problem, one can build the complete 
set of single-particle states consisting of the occupied states in the Fermi sea, 
h = (e/i, lh,jh, rnhi Th), and unoccupied states p = (ep, lp,jp, nip, Tp), where the 
indices Th{Tp) distinguish between neutrons and protons. Among the latter 
states we must distinguish between unoccupied states in the positive energy 
sector and those in the negative energy sector. We will use the notation p 
and p to specify them if needed. The occurrence of p states in RRPA is a 
consequence of the no-sea treatment of the mean field 1 10], 

The RRPA formalism is well-known and we need only to recall the main 
steps to make the notations clear. We denote by bp, 6]^ the creation operators 
for states p, h. The creation and annihilation operators of particle-hole states 
coupled to an angular momentum {JM) are: 
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where J = \J (2J + 1). In RPA, excited states \uJM) are created by the 
action of excitation operators Qt{JM) on the correlated RPA ground state 
\RPA). The excitation operators QI{JM) are hnear combinations of particle- 
hole operators: 

QliJM) = Y.Kt^ihiJM)+Y,Y;;;:'^cMJ-M), (2) 

ph ph 

where the X and Y amplitudes are determined by the RPA equations [11]. 

An important condition of the theory is that the correlated ground state 
must be a vacuum for all RPA excitations: 

Qu{JM) I RPA) = , for all vJM . (3) 

Then, it can be shown [12] that the RPA ground state can be built 
explicitly from the HF ground state \HF): 

I RPA) = iVoe^ |HF) , (4) 

where 

^ = ^ E E Mp,p,,,4,(JM)4,,,(J - M) , (5) 

JM php'h' 

and Nq is a normalization coefficient. The coefficients Mphp'h' are symmetric 
with respect to the ph and p'h' indices, Mphp'h' = ^p'h'ph- Note that, in all 
above equations the summations over p or p' run on both unoccupied positive 
energy (Fermi) p states and negative energy (Dirac) p states, as mentioned 
before. 

An important relation between the X, Y and AI coefficients can be de- 
rived if one uses the condition ^ that the RPA ground state is a vacuum for 
the Qy{JM) operators. This relation reads: 

Y.^php'h'4i) = -Y^l'\ (6) 

p'h' 

and it will be used in the next subsection. 

2.2 ONE-BODY DENSITY MATRIX OF RPA GROUND STATE 



Our main task in this work is to construct the one-body density matrix of 



the correlated ground state, p, 



'a/3 



RPA 

(0) 



that of the uncorrelated HF state, p L = ( HF 



RPA ) and to compare it with 



HF 
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We first observe that all quantities pph and php must vanish because 
operator products of the type and hpbh are just linear combinations of 
and Q, and neither (RPA |Qt| RPA) nor (RPA \Q\ RPA) can be different from 

zero. Thus, we are left to calculate the expectation values (RPA 



RPA^ 



and (RPA 



RPAy . The derivation was first given in Ref.|12j and we 
will not repeat it here. In spherical symmetry, the above expectation values 
are non zero only if the states p and p' , or h and /i', have identical quantum 
numbers except for their energies. The final form of the result makes use of 
Eq.®. We just give here the expressions necessary for calculating the RPA 
correlated density: 



RPA 



RPA\ = y-J%(fy;r^ 



vjh 



(7) 



RPA 



RFAj 



uJp 



(8) 



These results are formally identical with those of usual non-relativistic RPA. 
In the case of RRPA there is a major difference, however, as already pointed 
out in Subsection 2.1: the states p,p' can belong to the unoccupied Fermi 
states {p] or to the Dirac sea states {p}. 

In general, the corrections due to RPA correlations suffer from a double 
counting of the second order corrections, and therefore one should subtract out 
this redondant contribution. This double counting problem was analyzed by 
Ellis |13j who showed that it comes from the second order exchange diagrams 
of the RPA series. In the present case, however, we are dealing with a Hartree- 
type theory for the mean field, and the RPA series contains only diagrams of 
the direct type (summation of the bubble diagrams). Therefore, there is no 
need for evaluating double counting corrections in the present model. 

With the help of Eqs. ()7|8p and using the set of single-particle wave func- 
tions of RMF we can build the one-body density p{r) of the correlated ground 
state. The first term on the r.h.s. of Eq.([8]) gives rise to the uncorrelated RMF 
density p^'^\r), all other terms in YY correspond to corrections due to long 
range correlations of RPA type. 



3 CALCULATIONS 



In this work the calculations are performed with the effective Lagrangian DD- 
ME2 determined by Lalazissis et al.[9j for the RMF-RRPA approach. This 
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model gives a very good description of ground state properties at the (uncorre- 
lated) level, and also it accounts very well for the collective nuclear excitations 
in the framework of the RRPAj^. We limit ourselves to the case of spherical 
symmetry and we study the effects of long range correlations on the density 
distributions of the three magic nuclei ^^O, '''^Ca and ^^Ca. 

First, the self-consistent mean field is calculated by solving in coordinate 
space the RMF coupled equations for the nucleonic and mesonic fields. The 
detailed method is described for instance in Ref.[6j. This gives the complete 
set of single-particle states which will be necessary for the RRPA calculations. 
This complete set consists of the occupied states in the Fermi sea {h states), 
the unoccupied states above the Fermi sea (p states) and the empty states 
in the Dirac sea {p states). There are many p states which are unbound and 
therefore, we use a box boundary condition to discretize the continuum. The 
box radius R is chosen to be 12 fm for ^^O and 15 fm for the Ca nuclei. 

Next, the RRPA solutions are obtained by constructing the secular ma- 
trix ^ ^ ^ ) ™ particle-hole basis and diagonalizing it. The nec- 
essary formalism is explained in Ref.[7]. In this configuration space approach 
one has to adopt a truncation of the single-particle spectrum. In our case, we 
choose to keep all bound Dirac states whereas for the Fermi states we keep 
all states below Ecut= 120 MeV. Once the secular matrix is diagonalized we 
obtain the RRPA states u characterized by their amplitudes X^'^^ and Y^'^\ In 
the present work we include all the RPA states with natural parity vr = (— 1)"^ 
and angular momenta up to J = 8 in ^^O, J = 10 in ^''''^^Ca. 

One important property of RPA built on top of a self-consistent mean 
field is that the center-of-mass spurious state must be one of the J'^ = 1^ RPA 
states, and this spurious state must be at zero energy. In our calculations, we 
indeed obtain such a state and its RPA energy is very close (less than 500 
keV) to zero. It is well separated from other = 1~ states. We can thus 
clearly identify the spurious state and remove it from the states contributing 
to the expressions dZHHI). 

4 RESULTS AND DISCUSSION 

In this section we present and discuss the effects of RPA correlations on the 
densities of the doubly closed shell (or subshell) nuclei ^^O, ^''Ca and ^^Ca. 
According to Eqs.([7]l8]) the correlated densities p{r) can be decomposed into: 

p{r) = p'^°\r) + 6ppart.{r) + 5phoie{r) 

= p^^\r) + 6p{r) , (9) 
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where p^^\r) is the uncorrelated RMF density, and 



pp' vjh 

SpHoieir) = -Y.^Ur)^H{r)Y.PYX'^YX^^ ■ (10) 

hh' ujp 

Here, (ppir) and fhif) are the single-particle wave functions of the RMF. 
4.1 DENSITIES OF ^^O 

The respective contributions to Sppart. and 5phoie of the various multip clarities 
J'^ of the RPA states are shown in Figs|T] and [2j The calculations have been 
performed up to J'" = 8"^. The contributions to 6ppart. are generally posi- 
tive whereas those to dphoie are negative. The dominant contributions come 
from the multipoles 1~,2~^ and 3~. This is to be expected because the most 
collective RPA vibrations belong to these multipoles. Prom FigsH] and [2] it 
can be seen that the contributions decrease rapidly with increasing J and the 
convergence is reached at J'^ = 8"*". 
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Figure 1: Contributions to Sppart. of the J'^ multipoles from 0^ to 8^, in ^^O. 
The unlabelled curves correspond to = 6"'",7~,8+. 

Summing over the partial contributions we obtain the corrections 5ppart. , 
Sphoie and 6p which are displayed in FiglSl There are strong cancellations 
between 6ppart. and Sphoie-, and the net result is small. We note that the 
analytical structure of Eqs.dDlH]) insures that the integral of 5p over space is 
zero and hence, the total number of neutrons and protons is not affected by 
the RPA corrections. This is well verified in our numerical results. 
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Figure 2: Same as Fig.l, for 6phoie in ^^O. 
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Figure 3: The summed corrections Sppart., Sphoie and their total 6p in ^^O. 

In FigUJare shown the mass and charge densities, both in the correlated 
(RPA) and uncorrelated (RMF) cases. For the charge density the effects of 
the center-of-mass motion and of the finite proton size are included by the 
usual folding procedure [14J- For the mass density 5p{r) the correlations tend 
to fill up the hole in the central region while the surface becomes slightly 
steeper. The same tendency appears in the charge density to a lesser extent. 
In the inset of FigH] the charge density is shown in logarithmic scale in order 
to display the differences in the tails in the asymptotic region beyond 6 fm. It 
is seen that the correlations induce in the outer region a tail which decreases 
less rapidly than the uncorrelated density. To have a more global measure of 
the changes caused by the RRPA correlations, the root mean square (rms) 
radii of the densities are shown in Tabled) The changes on radii appear to be 
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modest, of the order of 0.7%-0.8%. 




r(fm) 



Figure 4: The uncorrelated (dashed hnes) and correlated (sohd hnes) mass 
and charge densities in ^^O. In the inset are shown the charge densities in log 
scale. 



Table 1: Mass and charge rms radii in "O and Ca, calculated with uncor- 
related RMF and correlated RPA densities. All units are in fm. 



Nucleus 










RMF 


2.594 


2.729 




RMF+RPA 


2.614 


2.747 


4UCa 


RMF 


3.346 


3.466 




RMF+RPA 


3.385 


3.501 



4.2 DENSITIES OF ^°Ca 

In FigsiSE] we present the respective contributions to 6ppart. and 6phoie of 
the various multipolarities J'^ of the RPA states, in the nucleus ^'^Ca. Here 
again, the dominant contributions come from the multipoles 1~,2~^ and 3~, 
and the convergence is reached at J'^ = 10^. It can be noted that the 1~ 
and 2+ contributions are localized near the nuclear center whereas the 3~ 
contributions extend over the nuclear volume. The summed corrections 5ppart. , 
Sphoie and Sp are displayed in FigJTl The net effect on 6p is still small because 
of the cancellations between Sppart. and 6phoie, but its shape is different from 
that of the ^^O case: the correction is negative from the center to about 3.5 
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Figure 5: Contributions to 5ppart. of the multipoles from 0^ to lO"*", in 
'^'^Ca. The unlabelled curves correspond to J'" = 6"^ , 7~ , 8"*" , 9^ , 10^. 

fm, and positive afterwards since the number of particles must be conserved. 
Thus, the change on rms radii becomes sHghtly larger and of the order of 
1%, as it can be seen in Table [H The mass density and charge distribution 
(corrected for center of mass and finite proton size effects as in the case of 
^^O) are shown in FiglHl 

4.3 DENSITIES OF ^^Ca 

For this nucleus with a sizable neutron excess we will examine separately the 
neutron and proton distributions and show that the long range correlations 
produce different effects on neutrons and on protons. In Figl9] are shown the 
partial contributions of various multipoles to 6p^p]^^ and & popart.- It can be 
seen that the different multipoles contribute in a similar way to the neutron 
and proton densities. Again, the dominant multipoles are those below J = 4 
and the convergence is reached before J = 10. The situation is slightly dif- 
ferent for Sp^j^^i^ and Sp^^J^^ as one can see in FigJlOl The most important 
difference between neutron and proton densities is seen in the contribution 
of the J'^ = 3~ multipole. This can be understood by looking at the ex- 
pression (fTO]l for Sphoie- In the case of neutrons the important configurations 
[(2sl/2, ld5/2, ld3/2)~^ (l/7/2)-^]3- are missing while they are allowed for pro- 
tons. This results in an important contribution from 3~ states to the proton 
density between and 2 fm, as compared to the neutron density. 

Summing up the contributions of various multipoles we obtain the density 
corrections Sp^^^^^l, Sp^j^^^J and their sum (Jp^^''^). The results are shown in 
FigfTTl for neutrons and protons separately. One can notice that the signs of 
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Figure 6: Same as Fig. 5, for Sphoie in ^'^Ca. 
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Figure 7: Same as Fig. 3, for ^"Ca. 

the corrections are opposite for neutrons and protons, in the central region up 
to 2 fm. The corrections become negative between 2.5 to 4 fm, and finahy they 
are positive in the outer part. In Fig {12] we show the mass and charge densities 
of '^^Ca calculated with the uncorrelated (RMF) and correlated (RMF+RPA) 
ground states. The general effect of correlations is to lower down slightly the 
density near the center and to shift the corresponding particles to the outer 
region. The exponent of the exponential tail is changed somewhat beyond 
8 fm. The values of rms radii for uncorrelated and correlated densities are 
summarized in Table O The changes in neutron and mass radii are relatively 
small (0.4% and 0.6%, respectively) whereas the changes of proton density 
radius (0.9%) and charge radius (0.8%) are more sizable. 

Finally, it is interesting to consider the respective roles of the Fermi and 
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Figure 8: Same as Fig.4, for ^"Ca. 

Table 2: Neutron, proton, mass and charge rms radii in ^^Ca, calculated with 
uncorrelated RMF and correlated RPA densities. All units are in fm. 



Nucleus 










4«Ca RMF 


3.576 


3.389 


3.499 


3.482 


RMF+RPA 


3.592 


3.420 


3.521 


3.511 



Dirac states in generating the corrections to the densities, in the case of RRPA 
correlations. In Figll3lwe show the separate contributions of the Fermi sector 
and Dirac sector to the density corrections Sppart., ^Phoie and their sum 6p. It 
can be seen that the contributions of the Fermi and Dirac sectors to Sppart. 
are both positive whereas they are both negative for Sphoif,. The magnitudes 
of Fermi and Dirac corrections are comparable. Thus, the Fermi and Dirac 
sectors play equally important roles in the density corrections. The fact that 
the total correction 6p remains small is due to the large cancellation between 
the positive Sppart. and the negative 6phoie- 

5 CONCLUSION 

In this work we have calculated the effects of RPA correlations on the ground 
state densities of some closed shell and closed subshell nuclei. This has been 
done for the relativistic mean field model and relativistic RPA theory, using a 
density-dependent effective Lagrangian with sigma, omega and rho couplings. 
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Figure 9: Contributions to 5ppart. (^^^ panel) and to Sp^part. (right panel) of 
the J'" multipoles from 0"^ to 10+, in ^'^Ca. The unlabelled curves correspond 
to J"" = 4+, 5", 6+, 7^,8+, 9^,10+. 



We have stressed the peculiarity of the relativistic approach as compared to 
the usual non-relativistic approach, with the important effects coming from 
the negative energy (Dirac) sector. The presence of the Dirac states is inherent 
to the relativistic approach because one needs the completeness of the single- 
particle basis. 

The long range correlation effects on densities are found to be limited 
although not negligible. It turns out that the main important multipolarities 
are the low J modes, which is to be expected because the important collective 
excitations correspond to these low multipoles. The partial contributions to 
densities decrease rapidly with increasing J, and it is sufficient to drop the 
contributions beyond J = 8 in ^^O and J = 10 in ^'^'^^Ca. 

We have used in this work the DD-ME2 effective Lagrangian which gives 
a good description of nuclear ground states in the RMF approach, and it 
is working well for collective excitations of electric type, i.e., transitions of 
natural parity vr = (— 1)'' described by RRPA without exchange terms[7j. 
However, it is not necessarily suitable for unnatural parity excitations, and this 
is the reason why we have not considered in this work the effects of magnetic 
transitions such as spin and spin-isospin modes. This must await future work 
where the effective Lagrangian is treated in Hartree+Fock approximation[15| 
and where the possibility to have pion- and p tensor-coupling would enable 
one to describe on the same footing the electric and magnetic transitions. 
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Figure 10: Same as Fig.9, for Sp^j^^^^ (left panel) and Sp^^^^^ (right panel) in 
48Ca. 
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Figure 13: The contributions of Fermi states (left panel) and Dirac states 
(right panel) to the mass densities 6ppart., Sphoie and 6p in ^^Ca. 
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